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The properties of the ferromagnetic frustrated spin-S one-dimensional Heisenberg model in the 
vicinity of the transition point from the ferromagnetic to the singlet ground state is studied using 
£SJ ■ the perturbation theory (PT) in small parameter characterizing the deviation from the transition 

point. The critical exponents defining the behavior of the ground state energy and spin correlation 
functions are determined using scaling estimates of infrared divergencies of the PT. It is shown that 
. the quantum fluctuations for s = 1/2 are sufficiently strong to change the classical critical exponents, 

while for spin systems with s > 1 the critical exponents remain classical. The dimerization in the 
singlet phase near the transition point is discussed. 
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I. INTRODUCTION 



The quantum spin chains with nearest-neighbor (NN) J\ and next-nearest-neighbor (NNN) J 2 interactions have 
been subject of numerous studies pj. The model with both antifcrromagnctic interactions Ji, J2 > (AF-AF model) 
is well studied @, IS 0, IE H II • The case of F-AF interactions (Ji < 0, J 2 > 0) is less studied. Though the latter 
model has been subject of many studies H, 0, 0, 0> E3 , the complete picture of the phases of this model as a 
function of the frustration parameter J%j ' J\ is unclear up to now. An additional motivation to study this model is 
• related to the fact that a class of recently synthesized compounds containing CuO chains with edge-sharing CuO^ 
units are described by the F-AF zigzag model [lj, [ijj, Ha li_Z| • The Cu — O — Cu angle in these compounds is close 
to 90° and usual antifcrromagnctic NN exchange between Cu ions is suppressed. This means that the sign of J\ can 
be negative, while the NNN exchange is antiferromagnetic. 
. It is well known that there is a critical value J2/J1 = —1/4, where the transition from the ferromagnetic ground 
state to the incommensurate singlet state occurs [l8l fl9| . The study of the character of this quantum transition is 
t-H , one of the interesting problem related to the F-AF model. In this paper we focus on the behavior of the model in 
the vicinity of the transition point. We hope that this analysis will be useful for the study of the properties of the 
edge-shared copper oxides where the frustration ratio is close to the critical point. In particular, for edge-shared 
cuprate LiiZrCuO\ the ratio is J 2 /Ji ~ —0.28 and for R02CU2M03O12 it is J2/J1 ~ —0.37 [20ll21j . 
The Hamiltonian of the F-AF model is 



N N 
H = -$^(S n ■ Sn+1 - S 2 ) + J^(S„ • S„ +2 - S 2 ) (1) 

where we put J% = — 1 and J 2 = J, s is a spin value and periodic boundary conditions are imposed. The constant 
shifts in Eq.QJ secure the energy of the ferromagnetic state to be zero. 

Unfortunately, this model is not solved exactly. As was noted above, the ground state of the model QJ is ferro- 
magnetic at < J < 1/4 and it becomes a singlet incommensurate state for J > 1/4. Though the transition point 
between these phases is J = 1/4 for any s 18], the spectra of the model with s = 1/2 and s > 1 in this point are 
different. For s = 1/2 the singlet ground state wave function in the transition point is known exactly |19L l22|. It is 
degenerate with the ferromagnetic state for any even N. For s > 1 the singlet ground state wave function is unknown. 
Finite-size calculations shows that at J = 1/4 the singlet state lies slightly higher than the ferromagnetic level and 
the energies of the singlet and the ferromagnetic states are equal in the limit N — > 00 only. 

In the vicinity of the transition point at0<7<Cl(7 = J— 1/4) the singlet ground state energy Eq behaves as 
Eq ~ ^ \ where (3 is a critical exponent. The classical approximation gives (3 = 2. The spin-wave theory as well as 
some other approximations |5l llOj do not change this critical exponent. Unfortunately, the exact diagonalization of 
finite chains shows a complicated irregular size dependence of the ground state energy, which makes the numerical 
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estimation of the critical exponent (3 impossible. In the paper |23| we conjectured that for s — 1/2 strong quantum 
fluctuations changes the critical exponent and (3 — 5/3. 

We note that the model Q with s > I has not been studied before and the critical exponent for these cases is 
unknown. In this paper we confirm our conjecture for s = 1/2 using scaling estimates of the perturbation theory (PT) 
in small parameter 7. We show also that /3 = 2 for s > 1, though the corresponding factor at j 2 is different from the 
classical value and it depends on s. 

One of the most important and open questions in the zigzag model {Q is the possibility of the spontaneous 
dimerization of the system in the singlet phase accompanying by a gap in the spectrum. This problem has been 
mostly studied in the limit of two weakly coupled AF s = 1/2 chains (J 1). The one-loop renormalization group 
analysis indicates [24I l25j that the gap is open. However, the existence of the gap has not been verified numerically 
pi). On the basis of a field theory consideration it has been proposed [2(| that a finite gap exists, but it is so tiny that 
it can not be observed numerically. On the opposite side of the singlet phase, J — > 1/4, there are no any reliable results 
about the dimerization and the gap. Strong nonmonotonic finite-size effects do not allow to study the dimerization 
numerically 

In order to study the problem of the spontaneous dimerization in the singlet phase of the model close to the 
transition point J = 1/4, we consider the generalization of the model Q by adding to the Hamiltonian H the 
perturbation in a form of dimerization term. Unfortunately, the used special version of the PT did not give us a 
rigorous answer about the spontaneous dimerization in the model However, it allowed us to estimate the critical 
exponent of the dimer order, in case if the spontaneous dimerization in the model Q exists. Besides, it allowed us to 
obtain the critical exponents of the ground state energy and the dimer order for the dimerized version of the model 

O- 

The paper is organized as follows. In Sec. II we present the scaling estimate of the critical exponent using the 
PT in 7 for the Hamiltonian (JTJ starting from the singlet ground state at J c = 1/4. In Sec. Ill we analyze the PT for 
the Hamiltonian which is transformed to new local axes forming a spiral structure. We establish the scaling behavior 
associated with 7 and with the pitch angle of the spiral. It is shown that the critical exponents for the ground state 
energy are different for the spins s = 1/2 and s > 1. In Sec. IV we study the problem of the spontaneous dimerization 
in the model Q. In Sec.V we summarize our results. 

II. SCALING ESTIMATE OF THE CRITICAL EXPONENT NEAR THE TRANSITION POINT J = 1/4 

We are interested in the behavior of the model in the vicinity of the transition point J c = 1/4. For this aim it 
is natural to develop the perturbation theory 

H = H + V-, 

H = -^(S„-S n+1 - s 2 ) + i]T(S ti .S„ +2 - s 2 ) 

V, = 7^(S„-S„ +2 - S 2 ) (2) 

with a small parameter 7= J-1/4<1 (7 > 0). 

At 7 > the ground state of the Hamiltonian H is a singlet. Since the perturbation V 1 conserves the total spin 
S 2 , the PT to the lowest singlet state |^o) of the Hamiltonian Hq involves only singlet excited states. The low-lying 
singlet excitations at the transition point have very small energies as shown in Figs.(l) and (2), where we present 
finite-size calculations of the energy gap between the two lowest singlets. These calculations show that the low-lying 
singlet excitations have different powers in N for the cases s = 1/2 and s > 1. As it will be shown below this fact 
leads to different critical exponents for spin systems with s = 1/2 and s > 1. 

The perturbation series for the singlet ground state energy can be written in a form: 

E ( 1 ) = (* \V 1 + V 1 V y + ...\* Q ) (3) 

Suppose that the low-lying excitations acting in the PT behave as 

E k -E ~ N- s (4) 

The higher orders of perturbation series contain more dangerous denominators, and, therefore, have higher powers 
of the infrared divergency. Therefore, we use scaling arguments to estimate the critical exponent for the ground-state 
energy. Below we will follow only powers of divergencies and omit all numerical factors. 
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Suppose that the matrix elements of the perturbation operator Vy between low-lying states involved into the PT 
at N — > oo behave as 

{^ l \V 1 \^,)^ 1 N 1 - d (5) 

with some exponent d. 

Looking after powers of infrared divergencies in all orders of the perturbation series the correction to the ground 
state energy takes a form: 



Eo(i) ~ (*o| |* ) E c ™ x ' m ~ ^ Nl ~ d ■ ( 6 ) 
where c 7ri are unknown constants and 



m— 



x ^ mvy\y ^ iNS+1 _ d (7) 

Ek — e 

is a scaling parameter, which absorbs the infrared divergencies. 

The scaling function f(x) at x —> is given by the first order correction. In the thermodynamic limit (x — > oo) the 
behavior of /(x) is generally unknown, but the natural condition £0(7) ^ -/V at iV — >■ 00 requires 

/(at) ~ x^+r^d (8) 

and, finally 

£0(7) 7V 7 4r^ (9) 

The perturbation series for the lowest excited state £1(7) has the same form as Eq.@. But a requirement of a 
finite mass gap (if any) m = £1(7) — £-0(7) ~ 0(1) leads to another critical exponent 

m ~ 7*+i-d (10) 

We note, that for the models in fixed points with a linear spectrum (S = 1), Eas.l(7|- (|lt)(l reduce to the well-known 
formulae [23] 

x = ~{N 2 ~ d 

£(7) ~ -JV-yS-H 

m — 7^ (11) 

where the exponent d represents the scaling dimension of the perturbation operator. 

However, the transition point J = 1/4 is not a fixed point. Finite size calculations for the gap between lowest 
singlet states give exponents (5 = 4 for s = 1/2 chain and 8 = 3 for s > 1 chain (see Figs.(l) and (2)). 

In order to determine the value of the exponent d, we notice that the singlet ground state of Ho has a spiral ordering 
at N — > 00 with a period of the spiral equal to N 

Oti-/ 

(*o|S„-S„ +i |*o) =s 2 cos— (12) 

For the case s = 1/2 this expression is an exact one ^3'E2j while for s > 1 we have observed the spin spiral 
structure in exact diagonalization of finite-size systems (the spin correlation function for spin s = 1 chain of size 
N = 20 is shown in Fig. 3). This means that the first order correction to the ground state energy in 7 is 

(*o|V 7 |* ) = -7^- (13) 
We assume that all matrix elements between any low-lying singlet states \^>i) and have the same iV-dependence 



mV-yl^)-^ (14) 



and, therefore, the exponent d = 2. 
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Thus, as follows from Eq.@ the critical exponents for the ground state energy are different for s = 1/2 and s > 1 

Eo(7) A^7 5/3 , s=l/2 

Eoif) ~ ~N-f 2 , s>l (15) 



For the case s > 1 the above scaling estimates reproduce the classical value for the critical exponent of the ground 
state energy. But in the special case s = 1/2 the quantum fluctuations are strong enough to change the critical 
exponent. In order to understand the nature of the difference between s = 1/2 and s > 1 systems and determine 
factors in Eas. (|15fl in the next Section we develop a special version of the PT. 



III. PERTURBATION THEORY FOR THE TRANSFORMED HAMILTONIAN 



Let us start from the classical picture of the ground state of the model Q . In the classical approximation the spins 
are vectors which form the spiral structure with a pitch angle tp between neighboring spins. The classical energy per 
site 

E c i(tp) = Ns 2 [1 - cos (p - J(l - cos(2^)] (16) 

is minimized by the angle 

(p cl = cos- 1 -^j (17) 

The classical ground state energy is 

2s 2 

£d(¥>ci) = -N— I 2 (18) 

Following this picture we transform local axes on n-th site by a rotation about the Y axis by tpn. This rotation 
transforms the original spin wave functions \4> n ) to a new basis depending on the angle if 

\*l>n, V ) = U v \lp n ) (19) 
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where 

is the rotation operator and C/J = J7_ v . For finite cyclic systems the pitch angle <p takes quantized values ip m — 
Under the unitary transformation U v the Hamiltonian H takes a form 

H v = U V HU- V =H + V V 

V v = (1- cosip) ]T [S:s: +1 + S z n S z n+1 ] 

-J(l - cos 2<p) [S X nS x n+2 + S z n S z n+2 ] 

- 5>in^(S*S* +1 - S z n S x n+1 ) - Jsm2 V (SZS z n+2 - S z n S x n+2 )] (21) 

Now let us choose some eigen state \ip n ) of the Hamiltonian H 

H\^ n ) = E n \^ n ) (22) 

The state \ip n ) is the eigen state of the Hamiltonian H, but not of H v . Therefore, if we develop and exactly calculate 
the perturbation theory in V v to this state we arrive to some eigen state \ip m ,ip) of the Hamiltonian H v 

H v \ip m ,ip) = E m (ip) \ip m ,tp) (23) 

corresponding to, generally speaking, another energy level E m (tp) ^ E n . Obviously, the unitary transformation U v 
does not change the spectrum. Therefore, the found energy level E m (ip) is also one of the eigen values of the original 
Hamiltonian H . Thus, taking different values of the pitch angle ip m — (m = 1 . . . N) and developing the PT in 
V v to some definite eigen state \ip n ) of the Hamiltonian H we obtain a set of N, generally different, levels E m (ip) of 
H. So, we do not need to fix the value of ip to its classical value in contrast to the spin-wave approximation. Instead, 
we are free to pick out the minimal energy from the set of the found N levels E m (<p). In the thermodynamic limit, 
when ip becomes continuous variable, this procedure means the minimization of the found energy E(tp) over (p. 
As a 'source' function \ip n ) of H it is natural to choose the ferromagnetic state with all spins pointing up 

|F) = |TT...T) (24) 
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FIG. 3: The spin correlation function (S« • S„+i) in the lowest singlet state at the transition point for spin s = 1 chain of size 
TV = 20. The spin spiral structure is obvious. 



This choice is equivalent to taking the function \FJ) — U V \F) as a probe ground state for the model QJ. The 
function \F ln ) has a spiral structure arising in the classical approximation. The expectation value of the total S 2 in 
this state is 

<^|S 2 |^Hy (25) 

This means that \F V ) is not a pure singlet state, but contains an admixture of states with S ^ 0. However, it is 
clear that the weights of states with 5^0 are negligible at N — ► oo and we can treat the state \F V ) as a singlet one. 

Since we are interested in the behavior of the model near the transition point J = 1/4, it is convenient to represent 
the Hamiltonian H v in the form 

H ip = H + V 1 + V v (26) 

with Ho and V 1 defined above in Eq.Q and to develop the perturbation theory to the ferromagnetic state in V — 
V v + Vj. So, there are two channels V v and V 1 in the perturbation theory characterized by two small parameters ip 
and 7. The ferromagnetic state \F) is the eigen state of the Hamiltonian Hq with the energy Eq — and also of the 
perturbation Vy, but not of V v . The obvious relation V 1 \F) =0 means that the perturbation series for the energy 
contains terms ~ <^ m 7™, but does not contain terms ~ 7™ without <p. 

At first sight it seems that as a result of the rotation Eci. (|2(J[) we obtain more complex Hamiltonian H v and more 
complicated perturbation theory with two channels. But the advantage of this method is to construct the perturbation 
theory in V = V v + Vy to the simple ferromagnetic state instead of the perturbation theory in V 1 to very complicated 
(and even unknown for s > 1) lowest singlet state of Hq, which was analyzed using scaling arguments and numerical 
calculations in the previous section. The fact that we separate the term V-, from H and treat it as the perturbation 
does not change our arguments about minimization of the found expression for energy E(ip,j) over ip. 

The ground state of the Hamiltonian Hq is manifold degenerate: all the ferromagnetic states \F$ Z } with different 
total S z = have zero energy. Therefore, at first we have to split this degeneration of the ground state with use 

of secular equation. It turns out that diagonal elements are proportional to (Fs z \ V \Fs z ) ~ N, while non-diagonal 
matrix elements are (F$ z | V \ Fs> z ) ~ 0(1). Therefore, in the thermodynamic limit we can neglect non-diagonal matrix 
elements and develop regular perturbation theory directly to the ferromagnetic state \F) with all spins pointing up. 

The first-order correction to the energy reproduces the leading terms of the classical result l|18f) : 

4 

£?W = (F\ V v \F) = -2Ns 2 j<p 2 + Ns 2 ^- (27) 
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The second-order correction to the energy 



relates to a two-magnon states, because operator V (actually V v ) have non-zero matrix elements in (|28JI only with 
the states \*5fk) containing two magnons with total quasi-momentum Q = and relative quasi-momentum k. Exact 
calculation of the two-magnon problem gives for the sum the following result 

E m = _ N 3sV (29) 

16(s + l) v ' 

This sum converges, because a dangerous denominator 



■s 



ea(fc) = E k - E = "-k A (30) 
for small k is compensated by the matrix elements in a numerator 

( * fc| ^ |F> = 4(JTTf 2 (31) 

As one can see, the two-magnon spectrum of Hq at Q = and 1 is simply twice an energy of one magnon 
E2{k) = 2ei(k) where 

£i(fc) = 2s(l - cosfc) - -(1 - cos(2fc)) (32) 

and E\(k) = sk^/A at small k. So, the low-lying states of Hq with small number of magnons have energies e m — 
mei(fc) ~ siV~ 4 , which leads to infrared divergencies in the next-order corrections to the energy. Similar to Eq.@ 
we sum them up using the scaling arguments. 

The PT for Eq. l|26|l contains two channels V 1 and V v , which are described by two independent scaling parameters. 
In order to determine these scaling parameters one should estimate large- N behavior of the matrix elements of the 
operators V 1 and V v between low-lying states and \*f?j), acting in the PT. Since the operators V 1 and V v create 
(annihilate) not more than two magnons, we look after only low-lying states with small number of magnons and 
energies 

e m ~ sN- A (33) 

We note that these states are very different from singlet states (with N/2 magnons) presented in Eq.© and this 
fact is crucial. 

The diagonal matrix elements for one-magnon states with small quasi-momentum k behave as 

(k\V 7 \k) = -Asjk 2 



(k\V v \k) = ^Vfc 2 ~^+W (34) 



(non-diagonal elements in the one-magnon sector are zero). 

For a small number of magnons m <C N we can treat them as almost independent, because the interactions between 
magnons gives only corrections of the order of magnon density p = m/N to the excitation energies and to the matrix 
elements. Therefore, large- TV behavior of the matrix elements (k ~ 1/-/V) are 



(tfilV^-) ~ i/sN- 2 (35) 



These formulae are validated by the exact solution of two-magnon problem. 

Now we are ready to identify the scaling parameters of the perturbations V 1 and V v . Similar to Eq.Q, they are 



m^m „ (36 ) 
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The scaling parameter ipN looks natural, because for the finite cyclic system the pitch angle tp is quantized as 

rm TV • 

The infrared divergencies are absorbed by these scaling parameters so that the divergent part of the perturbation 
series in both channels has a form 

oo 

E^ = (^V^) ]T c mn ((pN) 2m (yN 2 ) n (37) 

m,n=0 

with unknown constants c mn . 

In order to satisfy the thermodynamic relation E^ dlv ^ ~ iV, we rewrite Eq. (|37j) as 

£(*0 = Na<p * £ 9n (<pN) ( ) (38) 

where 

oo 

g n (<pN) = (<pN) 2n+2m - 3 (39) 

rn— 

are a set of (generally unknown) scaling functions. They should converge in the thermodynamic limit N — > oo to 
some constants 

a n = lim g n (ipN) (40) 

TV^oo 



which are Taylor coefficients of an unknown scaling function / 



Thus, collecting the converged first-order E^ and second-order E^ 1 corrections with the divergent part 

£(<M the 

energy takes the form 

s + 1 \ 



£ = -ZNs^v + Ns z Z- '— + Nsp> J f -3 (42) 



At 7 = the estimate of the energy relates to the spectrum of the Hamiltonian Hq. One can see that the energy 
for s > 1 is E ~ -/V(^ 4 , which for small <£> ~ ^ agrees with the numerical estimate ~ A" 3 (see Fig. 2). However, 
in the special case s = 1/2 the second term in Ea. H42|) vanishes and the energy becomes E ~ N<p 5 f (0), which for 
(/? ~ jj agrees again with the numerical estimate E ~ A^~ 4 (see Fig.l). From the positivity of the spectrum of the 
Hamiltonian Ho we conclude that / (0) > 0. 

Now we need to minimize the ground state energy over tp. As follows from Ea.l)42[l this procedure is different for 
s =1/2 and s > 1. For the case s = 1/2 

The comparison of powers in (p and 7 of two terms in Ea. H43l) shows that the minimum of E is reached at </?min ~ 7 1 / 3 . 
Therefore, / {^f— ^ — > / (0) at 7 — > and the expression for (^min takes a form 

VWfr) = (tU^) 1/3 (44) 



V5/(0) 

As was shown above / (0) > 0, which justifies Ea. H44|) . The corresponding minimal energy is 

E min = ^-^^7 5/3 (45) 

(2.5/(0)) 2/3 
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For the case s > 1 the minimum is defined by the first two terms in Ea. (|42|l and 



E min = SNs 2 ^!-^ (46) 
s — 1/2 

The last term in Ea. (|42[l gives the correction to the energy proportional to ~ Nsj 5 ^ 2 . 

Thus, we reproduce the critical exponents obtained in Sec. II. However, this special type of the PT allowed us to 
determine also the factor at 7 2 for the case s > 1, which at s — ► oo tends to the classical result Eq.JTSJ. 

According to Eas. l|44|l and l|46l) the pitch angle cp min has different behavior at 7 — * for s = 1/2 and s > 1. It does 
not coincide with its classical value l)17|l for any s, but it naturally tends to ip c \ at s — > 00. The found non-zero pitch 
angle ip m i n indicates the helical (spiral) structure of the ground state. Of course, this does not imply the helical long 
range order, which is destroyed by strong quantum fluctuations. Instead, this means an incommensurate behavior of 
the spin correlation function and the pitch angle </? m i n can be identified with the quasi-momentum q max at which the 
static structure factor takes its maximal value. 

IV. DIMERIZED ZIGZAG MODEL 

In order to study the problem of the spontaneous dimerization in the zigzag model Q close to transition point 
J = 1/4 we add to the Hamiltonian H the dimerization perturbation V a 

H d = H + V a 

V a = a]T(-l)™S„.S„ +1 (47) 
Then, the behavior of the ground state energy Eg(a, 7) of the model l|47(l gives us the dimerization order parameter: 

4 ( D _ irS „. S „ +1 ) = 4« (4 S, 

If £0(0:, 7) ~ —Napi^y) at a — > 0, then the singlet phase of the model is spontaneously dimerized and ^(7) is 
the corresponding order parameter. 

The classical approximation for the model 14711 shows that the spins form a double-spiral structure defined by two 
pitch angles tp and 9 so that the rotation angle about the Y axis on n-th site is 

(-1)" 

tp n = rup + (49) 
The expansion of the classical energy at (a, 7, (p, 9) -C 1 

E cl (a, 7, tp, 9) =N S 2 [^+ 9 ^- 2-fip 2 - atpej (50) 

is minimized by the angles 



tfd = V87 + 2a 2 

9 C \ = aip d (51) 

which gives the ground state energy at a, 7 C 1 

£ci(a, 7 ) =-^ 2 (4 7 + a 2 ) 2 (52) 

As follows from Ea. l|52|l E c \(a,j) vanishes on the line 

4 7 + a 2 = (53) 
which determines the transition line between the ferromagnetic and the singlet phases for the model l|47|l . 
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For the case s = 1/2 the exact singlet ground state on this line is known |2^,|2g. It has double-spiral long-range 
order 

(S, • Sj +n ) = i cos ip n 

where the angles ip n are defined by Ea. H49|) with pitch angle <p = ^ and small shift angle between spirals 9 = 

It is interesting that the classical relation 9 = cup (see Eq.lJSTJ) remains for the strong quantum s = 1/2 case on the 

transition line. The dimerization parameter on the transition line behaves as 

PtI = - = 7P a (54) 

Though p tI ^ the spontaneous dimerization is absent on the transition line in the thermodynamic limit. 
As follows from Eq. i|48[l the classical approximation yields the dimerization order for the model 147|l : 

p c i{a, 7) = sVci0 c i = 2s 2 a (4 7 + a 2 ) (55) 

Ea. l55|) shows that the dimerization vanishes on the transition line Ea . (|53|) . which accords with Ea.l|54[l. and it 
vanishes also at a = 0, which implies the absence of the spontaneous dimerization for the model {TJ. Since the 
classical approximation describes the limit s — > oo, one can expect that at least in the limit s — > oo the spontaneous 
dimerization in the model Q is absent. 

Following the classical picture we transform the local axes on n-th site by a rotation about the Y axis by angle 
(fin as written in Ea. l49|l . but not fixing <p and 9 to their classical values. Under this unitary transformation the 
Hamiltonian Hd l|47(l takes a form 

H Vt9 =H Q + V(a, 7 , <p, 9) (56) 

where the perturbation V(a, 7 , <p, 9) has a very cumbersome form and we do not present it here. 

Similar to the analysis done in Sec. Ill we develop PT in V(a,-f,ip,9) to the fully polarized state l|24() . The first 
order in V exactly reproduces the classical energy H5()(l 

=E cl (a,j,V,8) (57) 
The second order correction to the ground state energy gives 

E&> = -Ns 2 f V_ + V- a ^ 2 ) (58) 
U6(s+1) 2 J v ' 

As one can see the terms containing the angle 9 in the first order are exactly compensated by the contributions 
of the second order E^ 2 \ This result is rather unexpected. The classical approximation corresponds to the limit 
s — > oo and it would seem that the quantum effects will give relative corrections ~ s _1 to the energy. However, in 
this case the quantum corrections have the same order in s as the classical energy. 

The next-order corrections contain infrared divergencies and we treat them using the scaling arguments similar to 
that done in Sec. III. The analysis shows that the most divergent parts of the PT are accumulated in the following 
scaling parameters: 

x a ~ aN 
x v ~ ipN 

x 7 ~ 7A 2 (59) 

It turns out that the angle 9 is not accompanied by the infrared divergencies and, therefore, it does not form a 
scaling parameter. After the algebraic manipulations with the divergent series of the PT similar to Eq.^-lEl), the 
main contribution of the next-order corrections to the ground state energy at N — > 00 takes a form 



(60) 



where / (jjjs> is unknown scaling function of two scaling variables. The angle 9 does not contribute to the most 
divergent parts of the PT Eq.JHOJ), because it is not accompanied by the infrared divergencies. 
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Collecting the corrections E^ and E^ with the scaling part E^ tv ' we obtain the leading terms for the ground 
state energy: 

„2,,.. , j>j , AT 2<P 4 s - 1/2 , s , ( 1 a 



E = --Ns\A 1 + a )ip + JV 8 -r-__t_ + ATV/ S J (61) 

As follows from Ea. (|61|l the leading terms do not contain the angle 9. In fact, we have checked that the energy does 
not contain terms up to ~ # 4 . This result is not surprising. In general, the PT in V(a, 7, 92, 9) for the ferromagnetic 
state results in the energy E(a, 7, <p, 9) depending on 9. On the other hand, the spectra of the Hamiltonians Hd <|47[) 
and H Vt e ifSBT) coincide and the eigenvalues E n (a,^f) of both Hamiltonians do not depend on 9 and <p. Therefore, for 
any values of 9 and ip the PT leads to one of the determinate levels E n (a, 7) of the Hamiltonian H . The pitch angle 
ip is quantized as ip n = =jp- for finite N, and the PT with different ip n leads to generally different levels E n (a, 7). 
At the same time, in contrast to the pitch angle ip the angle 9 is a continuous variable even for finite N. Therefore, 
the continuity condition of the dependence E(a,j,ip,9) on 9 implies that the PT in V(a, 7, ip, 9) with any value of 
9 leads to the same energy level as at 9 = 0. In other words, the obtained in the PT ground state energy does not 
depend on 9. 

This fact is an argument for the absence of the spontaneous dimerization in the zigzag model (JTJ. Really, the PT 
with any value of 9 brings to the same state as it does at 6 = 0. But at 9 = and a — the PT in V(a, 7, 92, 9) 
reduces to the PT l|21|) in V v considered in Sec. III. There are no terms in the perturbation V v which break translational 
symmetry and can potentially lead to the dimer order. 

However, the above arguments do not prove the absence of the spontaneous dimerization in the model The 
rigorous method is to calculate the dimer order parameter directly from Ea. (|48|l . which we follow below. 

The minimization of the ground state energy over ip is performed in the same manner as was done in Sec. III. For 
the case s = 1/2 the second term in Ea. (|61|l disappears. The comparison of powers in ip of two terms in Eq.JfJJ) 

shows that the minimum of E is reached at tp m in ~ (47 + a 2 ) 1 / 3 . Therefore, we substitute / , ^7=^ / ^0, 

at 7 — > and the expression for <p 



(a, 1 ) = (4 1 + a 2 ) 1/3 g( V ) (62) 



-1/3 



where r\ = and g(rj) = [10/ (0, r])} 
The corresponding minimal energy is 



E min = -^(4 1 + a 2 f 3 g 2 ( V ) (63) 



For the case s > 1 the energy minimum is defined by the first two terms in Ea. H61|) and 

<Pmin = ^2 /_ + -J 2 V^7 + Q 2 

£min = -N S 2 ^ l {A 1 + a 2 ) 2 +N S {A 1 + a 2 f 2 g s {r 1 ) (64) 

where g s (77) = 2 ^zf/2 f s ^) ' ^ e l eacun § term in the ground state energy for s > 1 is determined by 

the regular parts of the PT, while the scaling part (the last term in Eq.JHU) gives only small correction to the energy. 

We see that the difference in the critical exponents for the cases s = 1/2 and s > 1 remains for more general 
dimerized model (|56|l as well. The pitch angle (pmin and the ground state energy _E m i n naturally vanish on the 
transition line Ea. (|53|l for both cases s = 1/2 and s > 1. 

The dimerization of the model l|56|l is defined as a derivative of the energy with respect to a (|48|l . As follows from 
Eas. l|63|l . (|64|) the dimerization at 7 = appears with critical exponents 

P\ 1= o - " 7/3 > « = 1/2 

Pl 7=0 ~ « 3 ' s > 1 (65) 

As for the model Q (<* = 0), the dimerization depends on the behavior of the scaling functions f s (77) at small 77. 
There are two possible scenarios. First, the expansion of f s (77) at 77 — > is f s (77) = a + brj^ with some constants a and 
b and fx > 1, so that f' s (0) = 0. In this case the dimer parameter is zero for the model QJ. Second, f s (77) = a + br) 
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at 77 — > and (0) = 6. For this case the translation symmetry of the zigzag model is spontaneously broken and 
dimer long-range order ^(7) appears as 

P\ a =o ~ 7 ?/6 , - = 1/2 

Pl Q =o ~ 7 2 , * > 1 (66) 

Here the critical exponent for the dimerization (|66(l for s > 1 comes not from the leading term in Eq. (|64|l , but from 
the scaling correction (the last term). 

Unfortunately, we do not have any information about the behavior of the scaling functions f s (rj). Therefore, we 
can only state that if the zigzag model Q is in the dimerized singlet phase at 7 > 0, then the critical exponents for 
the dimer LRO are given by Eqs.lJBBJl. 

V. SUMMARY 

We have studied the frustrated Heisenberg chain with the nearest ferromagnetic and the next-nearest neighbor 
antiferromagnetic exchange interactions. It was shown that the behavior of the model in the vicinity of the transition 
point between the ferromagnetic and the singlet phases depends on the value of the spin. For s = 1/2 the critical 
exponent characterizing the behavior of the energy is (3 = 5/3 in contrast to the 'classical' exponent = 2 for s > 1. 
This difference is a result of different finite-size dependencies of the spectrum at the transition point 7 = for the cases 
s = 1/2 and s > 1. The pitch angles characterizing the incommensurate behavior of the spin correlation functions are 
different for s = 1/2 and s > 1, too. In particular, the pitch angle <p of the spiral is proportional to 7 1 / 3 for s = 1/2 
and to 7 1 / 2 for s > 1. It means that the considered model with s = 1/2 is special and the quantum effects for this 
value of s are the most strong. 

One more intriguing question is related to the existence of the spontaneous dimerization in the singlet phase. In 
order to study this problem, we added to the Hamiltonian H the dimerization term and treated it as a perturbation. 
Unfortunately, the used special version of the PT did not give us a rigorous answer about the spontaneous dimerization 
in the singlet phase. Instead, under assumption of the existence of the spontaneous dimerization, the PT allowed us 
to estimate the critical exponent of the dimer order parameter. Besides, using the special version of the PT we obtain 
the critical exponents of the ground state energy and the dimer order for the dimerized version of the model ■ 
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